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Magnetars are slowly-rotating neutron stars with extremely 
strong magnetic fields (10? G [2]), episodically emitting 
~ 100 ms long X-ray bursts with energies of ~ 107974! erg. Rarely, 
they produce extremely bright, energetic giant flares that begin 
with a short (~ 0.2s), intense flash, followed by fainter, longer 
lasting emission modulated by the magnetar spin period (typi- 
cally 2 — 12s), thus confirming their origin [3] [4]. Over the last 
40 years, only three such flares have been observed in our local 
group [6] [4]; they all suffered from instrumental saturation 
due to their extreme intensity. It has been proposed that extra- 
galactic giant flares likely constitute a subset [7] [8] [9] of 
short gamma-ray bursts, noting that the sensitivity of current 
instrumentation prevents us from detecting the pulsating tail, 
while the initial bright flash is readily observable out to distances 
~ 10— 20 Mpc. Here, we report X- and gamma-ray observations 
of GRB 200415A, which exhibits a rapid onset, very fast time 
variability, flat spectra and significant sub-millisecond spectral 
evolution. These attributes match well with those expected for 
a giant flare from an extra-galactic magnetar [12], noting that 
GRB 2004154 is directionally associated with the galaxy NGC 
253 (—3.5 Mpc away) [13]. The detection of ~ 3 MeV photons 
provides definitive evidence for relativistic motion of the emit- 
ting plasma. The observed rapid spectral evolution can natu- 
rally be generated by radiation emanating from such rapidly- 


moving gas in a rotating magnetar. 


On April 15'^ 2020 at 08:48:05.563746 UTC, the Gamma-ray Burst Monitor 


(GBM) onboard the Fermi Gamma-Ray Space Telescope (Fermi) was triggered 
by an extremely bright, short and spectrally hard event, initially classified as a 
short Gamma-ray Burst (GRB), GRB 200415A [14], which was also detected by 
several other instruments [13] [15] [16]. An offline search using time-tagged event 
(TTE) data from the Burst Alert Telescope (BAT) onboard the Neil Gehrels 
Swift Observatory (Swift), obtained with the Gamma-ray Urgent Archiver for 
Novel Opportunities (GUANO) [T7] pipeline, also found the event. Using the 
light travel time of photons detected by the Inter-Planetary Network (IPN) of 
satellites, GRB 2004154 was triangulated to a 17 square arc-minute region cen- 
tered at RA and Dec. (J2000) of 11.88? (00h 47m 32s) and -25.263° (-25d 15' 
46"), respectively [I3]. The relatively small error box of the localization sig- 
nificantly overlaps with, and, therefore, is highly suggestive of, GRB 200415A 
originating from the Sculptor Galaxy (NGC 253), an active star-bursting inter- 
mediate spiral galaxy located ~3.5 Mpc away |18]. 

We use the BAT TTE data to determine the duration due to bandwidth 
saturation of the high-time resolution GBM TTE data (see Methods). We find 
the Too duration of GRB 2004154 (the time interval over which 5% to 95% of the 
total counts was accumulated [19]), to be 140.855 ms (1c). Correspondingly, 
the event Tso duration is 54.7*0 ms (1c). Further, our detailed temporal 


analysis of the event lightcurve shows that the rise time (10 to 9096) of the first 


pulse is Tii, = 77 + 23us (1o) (Fig.|1| panel (e)). 

We performed a timing analysis on the GBM light curve, searching for a 
rotational frequency in the range 0.02-50 Hz, finding no clear pulsation. We 
also searched the 40-4000 Hz window for Quasi-Periodic Oscillations (QPOs), 
possible signatures of seismic vibrations seen in the oscillating tails of confirmed 
Giant Flares (GFs); a candidate broad QPO was found at a frequency of v ~ 
180 Hz in the decaying tail of GRB 200415A with ~ 2.50 significance (see 


Methods). 

We performed time-integrated and time-resolved spectral analyses of the 
GBM data, focusing on the sub-ms structures in the lightcurve, shown in Figs. 
and The very high rate might cause the electronic signals of photons to 
overlap (pulse pile-up), causing their energies to be incorrectly measured, and 
spectral distortions. We evaluated this effect in the brightest interval (2 in 
Fig. |1) and determined that it was negligible, (see Methods). Among several 
spectral models used, we found that a power-law with an exponential cutoff 
(Comptonized) model, fit the data best; the Comptonized spectral parameters 
are presented in Table 1 (see Methods). The highest energy photons reliably 
associated with GRB 2004154, have energies of + 3 MeV (see Methods). Using 


time-resolved GBM spectral analysis with corrections from the BAT, we find a 


time-integrated isotropic equivalent energy output of Eis, = (1.51+0.021) x 1046 


erg (see Table 1). The peak isotropic luminosity is Lisomax = (1.53 + 0.13) x 


1055 erg s7}, while the total event luminosity is Liso = (1.07 +0.17) x 1047 erg 
s71, Our time-resolved spectral analysis shows remarkable sub-ms variations 
(Fig. Le) over à 10 ms interval, encompassing intervals 1, 2, 3, and part of 4. 
In Fig. 2d, we notice that the peak energy (Ep) reaches its highest value at the 
onset of interval 3, while it remains relatively constant throughout most of the 
event. 

The photon index (o), stays relatively constant at o ~ 0 during the event, 
which would be highly unusual for a short GRB. Figs. [2h,b show exponential 


decay trends in both energy flux F and Ep over interval 4, which is clearly 


discerned from the tail of GRB 200415A. The energy flux decay in Fig. [2]occurs 


on a timescale of rT = 45 + 3 ms. The E, is observed to decay on a longer 


timescale of r = 100 + 1 ms in the same figure. This exponential behaviour has 


been observed in other extra-galactic GF candidates [I3]. A distinctive F cx E> 


correlation was discovered (Fig. Dp: a signature of a relativistic wind. This 
unprecedented result is clearly observed in the GBM data for GRB 2004154, 
devoid of detector saturation effects. Such saturation effects likely precluded 
the ability to discern this trend cleanly from previous observations of galactic 
GFs from SGRs 1900+14 and 1806-20. 

Finally, we searched for radio emission associated with GRB 2004154, in 
four observations of the NGC 253 taken with the Karl G. Jansky Very Large 
Array (VLA) [20], 4.3 to 51.2 days after the event trigger. No significant variable 
or transient emission was identified. 

Previous studies postulated that about 1-20% of sGRBs could be extra- 
galactic GFs [7] 4 21] [22]. The sample of galactic GFs is very small and their 
properties are ill-determined due to instrumental effects from their extreme 
intensity. We therefore, first compare the GBM observations of GRB 200415A 


to the GBM observations of short GRBs (sGRBs) [24]. We find that the 64 ms 


peak photon flux Cpe ae = 73.7+2.1 photons cm ? s^!) of GRB 200415A lies 
at the 97.5" percentile of the sGRB distribution, the peak energy (Eg^t^les = 


9984-45 keV) at the 79t” percentile, and the photon index (a*?t*2le$ = 0.39--0.09) 


at the 88.5} percentile. It is similarly near the edge of the a distribution for 
the burst population detected with the Burst and Transient Source Experiment 
(BATSE) on the Compton Gamma-ray Observatory [25]. Consequently, we find 
the flat, hard spectral slope, high Ey and peak flux during the brightest 64 
ms of GRB 200415A (shown in Figs. [1| and [2) to be unusual for sGRBs, thus 
better explained as the initial spike of a magnetar GF from NGC 253. This 
interpretation is further motivated by similarities of the properties of this event 
to previously proposed extra-galactic GF candidates [10] 23]. A rapid rise time 
is characteristic of a GF onset. Compared to the rise times reported in the 


GRB catalogs of GBM and BATSE, this rise time is considerably shorter than 


for any event in these samples, and is shorter than extreme examples of short 
variations, such as ~ 100 us for GRB 910711 [26] and 2.8 ms for GRB 090228 
E. 

Unfortunately, we could not detect the magnetar period-modulated tail with 
putatively E ~ 10% erg. This ‘smoking gun’ evidence for a GF, is observed over 
several hundreds of seconds in all three confirmed magnetar GFs, but is absent 
in GRB 2004151 as it is likely below the detection threshold for GBM given 
its distance to NGC 253. For other extra-galactic GF candidates this feature is 
similarly undetected [13]. 

'The standard picture for the origin of the GFs is the release of energy trig- 
gered by fracturing the crust of the neutron star[1] by large sub-surface magnetic 
fields, depositing hot plasma into the inner magnetosphere. Using a GF interpre- 
tation, the GBM observations indicate that the MeV-band emission must come 
from a relativistic outflow that is initially, highly ‘optically thick’ (opaque). 
The enormous isotropic equivalent luminosity of Liso > 10^" erg s~! is orders of 
magnitude larger than the fiducial Eddington luminosity limit, Lgaa ~ 10% erg 
5-1, for a neutron star of solar mass, Mo [28]. This limit defines when radia- 
tion pressure associated with electron (Compton) scattering overwhelms gravity 
and pushes hydrogenic gas away from the surface to high altitudes. For GRB 
2004154, we thus expect a relativistic wind [1] with bulk Lorentz factor T > 1 
(i.e., speed c(1— 1/T?)!/2) to be present, putatively over the magnetic poles. At 
high altitudes, the radiation pressure abates and the wind “coasts” with con- 
stant I. Transparency of the wind to quantum electrodynamical (QED) pair 
production yy — ete" of electrons (e^) and positrons (e+) by 3 MeV photons, 
detected by GBM, unambiguously implies that T > 6 (see Methods). This is 
a more stringent bound than was possible for the < 1 MeV photons seen in 


the initial spikes of galactic GFs [6] [29]. This GBM limit is consistent with the 


stronger constraints due to the detection of GeV photons by Fermi-LAT [15]. 
High Lorentz factors (T Z 100) are predicted from dynamical models [30] of 
thermal fireballs with high peak energies, usually applied in GRB emission mod- 
eling. The high opacity to QED magnetic pair creation y — e*e^ in the inner 


magnetosphere of magnetars [31] indicates that this wind is likely dominated by 


e pairs, with limited baryonic content. The dense wind relativistically boosts 
its embedded radiation to higher frequencies via the Doppler effect and “beams” 
or collimates it into a radiation emission “cone” of opening angle Ocon ~ 1/T 
radians. The observed correlation between the energy flux and the peak energy, 
F x E? (Fig. |2) can be readily explained by relativistic Doppler boosting (see 
Methods). 

The GBM spectrum is very flat, ie. has a relatively high value for its 
power-law spectral index a. This is expected for a wind that for the most part 
is highly opaque to electron scattering, a so-called (ballistic) Compton cloud. 
The huge Liso indicates electron number densities of ne = 102 cm~? when the 
plasma is close to the stellar surface, for which the electron-photon Compton 
scattering mean free path is A, ~ 10-7? cm or shorter. Opacity to scattering 
persists out to altitudes of R ~ 10? — 10!? cm and higher. The very flat indices, 
a œ% —0.2 — 0.3 are similar to those identified in Lin et al.[32] for normal 
SGR J1550-5418 bursts, wherein it was concluded that the high ne densities 
must seed rampant Comptonization that yields spectra beginning to approach 
a modified blackbody (Wien) form [28]. The same situation is anticipated for 
GRB 2004154, with its markedly higher deduced n, values. Yet, the radiation 
is not truly thermal: the photon energy flux at the source as inferred from the 
GBM-detected flux is substantially inferior to that appropriate for a Planck 
distribution from a relativistic wind that yields the same Ep in the observer 


frame (see Methods). This strongly contrasts the picture of GRB fireballs, which 


possess fluxes and Liso many orders of magnitude larger, yet with similar Ep 
values. The broad, flat spectrum of the GF may in fact be a superposition 
of Comptonized Wien-like spectra from different altitudes spanning a range of 
effective temperatures in the adiabatically cooling wind. Note that such flat a 
indices are inconsistent with synchrotron emission scenarios commonly invoked 
for GRB spectral interpretation [33]. 

The spectral evolutionary sequence displayed in Fig. |1| suggests that the 
initial hardening might reflect relativistically-boosted emission coincident with 
later stages of the wind acceleration epoch, followed by some cooling during the 
coasting phase of the laterally-expanding wind. Alternatively, it could be a con- 
volution of the influences of stellar rotation and radiation beaming (collimation) 
into a “cone” subtending an angle of Ocon ~ 1/T' . As the Doppler-boosted cone 
sweeps across our line of sight, the intensity would thus increase, peak and then 
decline, accompanied by a spectral hardening and subsequent softening; this is 
commensurate with the evolution displayed in Fig. Such a transient “rela- 
tivistic lighthouse beaming” effect would generate the decay time 7 & 45 ms 
of the tail in Fig. [2] of a P = 8s rotation period for a IT ~ 30 bulk flow (see 


Methods), noting the intrinsic coupling 7 ~ P/(2aT) for the rotating beam. 
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Figure 1: Temporal and Spectral Variability of GRB 200415A. (Left:) 
Lightcurves with 0.2 ms resolution for a Fermi GBM BGO detector (a), Nal 
detector (b) and Swift-BAT (c). The BAT lightcurve was shifted by 5.7 us to 
account for the light-travel time between the spacecrafts. Panel (d) shows the 
spectra for the four intervals. The shaded area indicates 1o confidence regions. 
Using the BAT TTE data, we identify that the GBM TTE bandwidth (see 
Methods) was exceeded from -2.4 ms to -0.8 ms (horizontal black lines in panels 
a,b), resulting in a 47.3% loss of flux in interval (2). There is also a ~1.8 ms 
data gap from 4.6 to 6.4 ms caused by a CSPEC packet blocking the GBM TTE 
data, resulting in a 3.47% loss in interval (4). Panel (e) shows the first pulse 
with high temporal resolution (12 us) and the fitted pulse profile. 
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Figure 2: Flux and Spectral Evolution of GRB 2004154. Spectral vari- 
ability of GRB 200415A. Panels a, b and c show the energy flux or F (a), peak 
energy or E, (b) and a (c) evolution over the Too duration (~141 ms). t is 
the time in seconds. The blue lines are exponential fits to F and E, over the 
tail emission of GRB 200415A (5 ms after zero time). This results in a decay 
timescale of 45 + 3 ms for F. The decay time for E, is 100 + 1 ms. The tem- 
poral binning is 8 ms for these panels. Panels d and e show how the E, and a 
trend on sub-ms timescales (At = 250 ys), over the interval shown between the 
dashed grey lines in Panels a through c. Panel f, shows the relationship between 
F and E, over parts of the burst not affected by data saturation, using 8 ms 
temporal resolution. The blue line is a power-law fit to this data (exclusively 
after the main peak), which shows F œ EŻ. The spectral range used for these 
measurements was 8 keV to 10 MeV. The coloring scheme follows that in Fig. 
The zero time reflects the time of the GBM trigger. All fit errors are at the lo 
confidence level. For more information on their derivation, see Methods. 


Methods 


GBM Observations and Data Processing Fermi-GBM consists of 12 thallium- 
doped Sodium Iodide (Nal) detectors and two Bismuth Germanate (BGO) 
detectors. The uncollimated Nal scintillator detectors are clustered into four 
groups of three detectors at each corner of the spacecraft, arranged such that 


any cosmic source above the Earth's horizon will illuminate one cluster. The 
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combined effective spectral range of both the Nal and BGO detectors are + 8 
keV to 40 MeV. GBM has several data types, which are either produced by a 
trigger or continuously. Triggered data types include CTIME data (binned to 
64 ms resolution with 8 energy channels), and CSPEC data (binned to 1024 
ms with 128 energy channels). The GBM TTE triggered data are the tagging 
of individual counts in the detectors within 5 minutes of the trigger time and 
have a relative timing resolution of ~ 2 us and 128 energy channels. The energy 
channels for all of these data types are pseudo-logarithmically spaced. More 
information on these and other data types, as well as other details pertaining 
to the instrument, can be found in Meegan et al.[34]. 

A bandwidth limit between GBM and the Fermi spacecraft restricts the 
time-tagged GBM event rate to ~ 375,000 events per second, summed over all 
detectors and energies [35]. When the GBM TTE rate exceeds this value, TTE 
counts are lost in a manner that is not biased by detector or energy (see Fig. 
1), consequently spectral analysis to obtain spectral shape parameters is still 
possible, but the flux normalization will underestimated. During the saturated 
interval, sub-intervals may be weighted differently, but the data is still in correct 
order. Additionally the transmission of another GBM datatype to the spacecraft 
can block the T'TE data, causing a gap. Other GBM data types are not effected 
by this bandwidth limit but have inadequate temporal resolution for this event. 

Due to high count rates, the electronic signals of the counts can overlap 
(pulse pile-up), causing incorrect energy measurements and spectral distortions. 
To check for this in the GBM data for GRB 2004154, we used an analytical 
method [40] that was verified by high-rate measurements with a GBM detec- 
tor [39]. This method was applied to the spectral model and rate of interval 2, 
for the detector with the highest flux (Nal 1). The spectral shape of the detector 


counts are only slightly modified by pulse pile-up, with a change in slope below 
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30 keV and a loss of higher energy counts starting at 200 keV, reaching a 5% 
loss above 400 keV. We note that in joint fits, the lower-rate BGO data help 
constrain the results. Pulse pile-up is otherwise not included in our analyses. 

The localization accuracy of an all-sky monitor like GBM is limited, espe- 
cially for short events. While an initial position for GRB 200415A was promptly 
determined using GBM [86] [37] using several methods [38] 41], each position 
had several degrees of error. We thus used the IPN localization[13] to generate 
Detector Response Matrices (DRMs). 

Each GBM detector was checked to see whether their viewing angles were 
within 60? of their respective on-axis position to the source, and whether they 
were blocked by parts of the spacecraft. From this analysis, we identified Nal 
detectors 0, 1, 2, 3 and 5 and the BGO 0 detector as having satisfied these 
criteria. Spectral parameters were determined by selecting the model that best 
fit the data, by looking at variations in the Cash-Statistic (C-Stat.) per degree 
of freedom (dof) [42]. A simple power-law, a power-law with an exponential 
cut-off (Comptonized), a black-body, a Band function and combinations of 
these models were used. The values in Table 1 and in Figs. [1| and [2] are all 
spectral parameters derived form the Comptonized model fit to the data, which 
was found to be the best. The estimated covariance of the optimal values to the 
fits of the spectral parameters in Fig. [2h.b,f, use the minimization of the sum 
of the squared residuals. The diagonals provide the variance of the estimated 
parameter. One standard deviation error on the parameters was done by taking 


the square-root of the diagonals from the covariance matrix. The full fits to the 


data in Fig. 2h.b.f, are F = (9.39-E0.39) x 1075 . e(-0.022:0.000*. E. = (1299 + 


82) + e(-0-01040.001) and F = (5.8440.58) x 10-5 - (E,/1000 keV )2-0450:37, 


respectively. F is measured in erg cm? s~!, t in seconds and E, in keV. 


All of which include 10 errors, as shown in the Fig. This results in 
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decay times (with lo errors) of 45 + 3 ms and 100 + 1 ms for F and Ep, 
respectively. All subsequent reported times are measured relative to the trigger 
time at 08:48:05.563746 UTC on April 15^^ 2020. 

As a precursor was observed 142 s before the initial spike of the galactic GF 
from SGR 1806-20 [29], we searched for precursor emission in a 2400 s interval 
preceding GRB 200415A using the GBM targeted search [44]: this search 


returned no candidates. 


Swift-B AT Observations and Data Processing The BAT is a coded-mask 
imaging instrument with a 2 steradian field of view (FoV) in the 15-350 keV 
energy band [45], but is capable of recording photons with energies up to ~ 500 
keV. For GRB 2004154, we used a new capability to provide data from BAT: 
the Gamma-ray Urgent Archiver for Novel Opportunities (GUANO) [I7], which 
allows the time-tagged BAT event data to be saved on demand, based on triggers 
from external instruments. At the time of the burst, BAT was oriented such 
that the source illuminated the detector at an angle of 48.51? from the center 
of the FoV. Although the location was entirely outside of the coded FoV (by a 
few degrees), the source was sufficiently bright to penetrate the graded-Z shield 
on the sides of the instrument and deposit significant flux onto the detector 
array. If the burst had been within the FoV, it would have resulted in a prompt 
autonomous trigger and most likely, would have saturated the instrument. 
Using an instrument and spacecraft mass model (Swift Mass Model: SwiMM 
[46]) in GEANT4 and Monte Carlo methods, we generated custom DRMs for 
the orientation of the source, thus allowing the use of XSPEC to fit spectra 
and derive fluxes from the count rates. The effective area response of the BAT in 
this orientation is both substantially reduced and altered with respect to sources 
originating from within the FoV, due to the shielding. In this orientation, the 


BAT effective area peaks at 200 keV with an effective area of ~ 400 cm?, 
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quickly dropping to « 50 cm? at around 50 keV. 

The BAT data were also used for timing analysis, such as the Too and T5 
durations, shown in ED Fig. 1. The errors for this measurement were estimated 
by taking a standard deviation of the median fluence value after the cumulative 


counts curve flattened. 


Rise time We perform an analysis of the first pulse using the GBM TTE data 
from the Nal and BGO detectors (interval (1), Fig. [i{e)), that starts at —4.4 
ms and is unaffected by data saturation. We fit a pulse function introduced by 
Norris et al. [48] to model the properties of GRB pulses. The count rate as a 
function of time is described by I(t) x exp(—((Zpeak — t)/Ti)”) before the peak 


of the lightcurve and I(t) x exp(—((t — Tpeak)/T2)”) after the peak. We find 


Ti = 32.8 + 9.945. The rise time, defined as the time elapsed between the 10% 


and 90% of the peak yields Tiise = 77 + 23us at the lo confidence level. 


Energetics and Time-resolved Spectra We divided the lightcurve into 4 
intervals (see Fig. |l| and Table [1]. 'The differential photon number spectrum 
(dN/dE) in all four intervals is best described by a power law with an exponen- 
tial cutoff. The cutoff is parametrized as the peak energy of the vF, spectrum, 
Ey. The first interval (1) contains the first pulse, present in both GBM and 
BAT. It is a ‘clean’ pulse, not affected by saturation and there is no overlap with 
other pulses. Interval (2) is the brightest part of the lightcurve, affected by TTE 
saturation in GBM. Interval (3) is the hardest with peak energy Ey — 1.9 MeV. 
The fourth interval has a featureless decay out to 136.4 ms and contains most of 
the fluence. The spectral fits for intervals (1), (3) and (4) are illustrated in ED 
Fig. 2. The comparison of the energetics to SGRBs were over 64 ms, a timescale 


generally reported with sGRB properties in mission catalogues. 


Highest Energy Photons ED Fig. 3 shows the individual TTE counts in 


GBM BGO detector 0. The GBM detectors register energy deposits and are 


19 


unable to identify particle type, nor can they determine photon arrival direction. 
Therefore, it is impossible to determine with certainty whether any particular 
TTE is from GRB 2004151, another gamma-ray source or is due to other back- 
ground. Instead, we can determine whether a rate increase is statistically signif- 
icant and therefore associated with GRB 200415A. We use a Bayesian method 
applicable to Poisson rate data, for the classic on-source / off-source method of 
source detection [49]. The method uses data from two time intervals, an off- 
source or background interval and an on-source interval, to test two hypotheses: 
1) that all of the T'TE are due to background, and 2) that there are excess TTE 
above background in the on-source interval due to the source. The method 
needs a prior expectation for the source count rate, which we obtain from the 
spectral fit in the on-source interval. For the 2.5 to 3.5 MeV energy range, there 
are 100 TTE in a 0.99 s background interval and 9 TTE in a 6.4 ms on-source 
window (blue box in ED Fig. 3 The calculated probability for a source signal 
(GRB 200415A) above the background is 0.9999997. This energy range is well 
above the threshold for yy pair creation; the two TTE with the highest energies 
are 3.0 and 3.1 MeV. 

We also consider a higher energy range, 3.5 to 10 MeV. In this energy range 
there are 88 off-source TTE and three on-source TTE (red box in ED Fig. 3) 
with energies of 4.0, 6.7 and 8.8 MeV. The probability that the three on-source 
TTE are an excess rate that should be attributed to GRB 2004154 is 0.966. 
Consequently, we do not consider these three TTE sufficiently significant to use 
for our Lorentz factor analysis. 

The definitive detection of Emax ~ 3MeV photons for GRB 200415A by 
GBM is at an energy well above the pair-production threshold. Dedicated 
gamma-ray instrument observations of the main pulse of galactic GFs have typ- 


ically been saturated, and likely suffer from spectral distortions (pulse pileup). 
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Thin particle detectors have reported usable spectra (e.g., the WIND SST sili- 
con detectors measured a spectrum for the December 27, 2004 GF main pulse 
of SGR 1806-20 extending to 1 MeV [29]). Gamma-ray detector observations of 
extra-galactic GFs will likely provide the best results on their spectra, obtaining 


excellent statistics without spectral distortions due to extreme fluxes. 


The Environment of the MeV Emission Region The appearance of emis- 
sion above the two-photon pair creation (yy — ete” ) threshold of mec? = 511 
keV in GRB 200415A can be used to provide a lower bound to the wind bulk 
Lorentz factor [ relative to the magnetar. This is a traditional practice in GRB 
studies [50] [51] [52], using the bulk motion to reduce the pair opacity through rel- 
ativistic collimation of the radiation field within an angle Ocon « 1. The spec- 
trum of this giant flare with its emission limited to energies E < Emax ~ 3 MeV 
demands an individualized calculation of pair opacity. The most conservative 
estimate corresponds to all observed GBM photons being below the 511 keV 
threshold in the comoving frame of the plasma/photon gas. Then one sim- 
ply has T > Emax/ mec? ~ 6. Somewhat higher values can be obtained with 
additional assumptions of source spectral extension beyond Emax, albeit unde- 
tectable due to the radiation background. The pair opacity constraint on T is 
more restrictive than is possible for Emax ~ 1 MeV values appropriate for the 
SGR 19004-14 and SGR 1806-20 giant flares. Expectations from isotropic fire- 
ball dynamics modeling and especially Fermi-LAT observations suggest 
much larger values for the bulk T. 

'The GBM emission is non-thermal. For a distance d — 3.5 Mpc to the Sculp- 
tor galaxy, the detected photon energy flux can be mapped over to an energy 
flux F(R) through a surface located at a distance R from the magnetar's center. 
The entries in Table[1]for interval (3) indicate that F(R) ~ 2 x 109? (105cm/ R)? 


erg em ?s- 1. If this radiation were to be thermal, then this flux should be com- 
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parable to the Doppler-boosted Planck spectrum that generates the energy peak 
value Ey ~ 1856 keV for this interval. Suppose the Planck spectrum in the wind 
frame has some temperature Tw. For a wind moving with speed $c and bulk 
Lorentz factor  — (1 — 6?)~1/? > 1, this yields a temperature Tog = dpTw ~ 
E,/3k in the observer's frame. Here dp = [I'(1 — 8 cos655,)] ! is the Doppler 
factor of the wind for an observer viewing it at an angle Oops to its velocity. Gen- 
erally, the most intense emission arises when 4%, < 1/T is small and the wind 
is viewed head-on, corresponding to óp ~ DP > 1. The energy density of the 
thermal radiation [28] in the wind frame is Uj, = 72/15 (m,c?/M3) O1 = aT’, 
for dimensionless temperature Oy = kTy,/mec?, where A, = h/mec is the re- 
duced Compton wavelength of an electron, and a = 7.56 x 107 !?erg cm? K~+ 
is the radiation constant[28]. The corresponding pressure is Pin = U,5/3, and 
its enthalpy is Win = Um + Pin = 4Uin/3. Boosting to the star frame, the 
energy flux of the photons through a surface locally perpendicular to this boost 
is JAn (0p) & 623Winc = 4n7/45 (mec? / M3) 62, 01, so that Fin(dp) x 62. We 
choose bulk dp ~ I ~ 100, as deduced from simple spherical dynamics argu- 
ments [15]. Setting Ow — kTe¢/(dpmec?) quickly yields a thermal energy flux 
JA (6p) ~ 8 x 1091 (102/605)? erg em ?s-!. This is clearly much larger than the 
aforementioned flux inferred from observations, and much more so assuming 
that the radiation emanates from altitudes R > 10? — 10!?cm where the plasma 
becomes transparent to electron scattering. Adopting a smaller Doppler factor, 
dp ~ 10, accommodating more closely the T ~ 6 bound obtained from the pair 
transparency considerations, increases the thermal flux and renders the dispar- 
ity with the observed flux more extreme. Thus, the thermalization in the wind 
is incomplete, consistent with the inherently non-Planckian form of the spectra 
at all times. Then, the radiation pressure exerted on the plasma is inferior to 


that invoked by applying traditional spherical GRB thermal fireball models [53], 
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possibly by factors of ~ 10 — 104. 

The GBM emission is subjected to prolific Compton/Thomson scattering 
by electrons and positrons, which shape the non-thermal spectrum. The ob- 
server frame pair density ne of the wind couples to the radiation density ny 
through a multi-layered interplay between expansion dynamics, geometry of flar- 
ing magnetic field lines that guide the outflow at altitudes below 108 — 10? cm, 
radiative transfer in a strongly-magnetized plasma, and magnetic opacity and 
pair equilibria in the wind frame. The complexities of these are beyond the 
scope of this paper. Yet, for a mean electron energy l'm.c? in the observer 
frame, one can connect n, to the emitted energy flux via a simple estimate 
en, (Pm,c?) c ~ F(R) for some unknown radiative efficiency e that is anticipated 
to not be vastly different from unity. Using Table 1 for interval (3) indicates that 
F(R) ~ 2x10?*(1065cm/ R)? erg cm-?s-!. Assuming e ~ 1, one can discern that 
ne ~ 8 x 107? cm"? at the stellar surface (where T ~ 1), and the optical depth 
Tr = ne04 HR to Thomson scattering (of cross section or = 6.65 x 10-75 cm?) 
is probably ~ 5 x 10! there. Tr doesn't drop below unity until altitudes of 
R > 10!°cm, since n, x R^?, because magnetic field line colatitude 0 and ra- 
dius R satisfy sin? 0  R in flaring dipolar field geometry. Magnetic reduction 
of the scattering opacity can lower this transparency radius by a decade or 
so, and concomitantly decrease the radiation pressure and lower the ultimate 
bulk Lorentz factor somewhat. 

À key GBM result is the exponentially decaying tail of the flux in Interval 
4 (3-136.4 ms), after the initial spike (see Fig. ). During this decay, which is 
of a timescale, T ~ 45ms, there is a distinctive F œ E? correlation, depicted in 
Fig. This is the hallmark of a relativistic wind. The intrinsic couplings of 
observer frame flux F œ 02 and peak energy Ep x dp identified above, are true 


for a relativistic boost from the wind frame regardless of whether the photons are 
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thermalized in the wind or not. Combined, they naturally generate the observed 
F x ES correlation at an approximately fixed emission radius. A distribution of 
Doppler factors dp is anticipated when sampling flare evolution as the observer 
angle fobs changes significantly as the magnetar rotates, precipitating spread in 
the E, and F distributions. Flaring of the magnetic field lines when sampling 
different emission radii may modestly broaden these distributions. 

The extremely short rise time of Tii. ~ 80s corresponds to a physical 
scale of cT, ~ 2.4 x 106cm. The tail decay time is a factor of ~ 550 longer. 
These can potentially constrain emission region sizes that are vastly smaller 
than those inferred for GRBs from their variability times. Yet this connection 
is impacted by the physical rotation of the magnetar. The tail decay timescale 
of T = 45 ms corresponds to a stellar rotation through an angle of A0 = 2? for a 
magnetar of typical period (P — 8s). This type of duration is naturally expected 
for relativistic beaming of radiation from outflows with T ~ 1/@con ~ 30 for 
Oco ~ 2°. Even if the physical angular extent of the emission region exceeds the 
relativistic collimation angle O. oy, Doppler boosting dictates that the dominant 
signal will arise when the instantaneous observing angle to the wind velocity is of 
the order of 055, < Ocon. Thus, a picture emerges that what we might be seeing 
in the GBM, Swift-BAT and KONUS-Wind signals is produced through stellar 
rotation of the Lorentz of Doppler-boosted emission in and out of our view. 
Accordingly, temporal scales At ~ P/(2xT') may be signatures of the relativistic 
beaming structure associated with the GF wind convolved with the magnetar's 
rotation, what can be termed a relativistic lighthouse sweeping effect. The 
greater Doppler boosting near the Oops < Ocon “core” of the Lorentz cone would 
yield harder spectra near the peak of intensity and softer in its ingress and 


egress, as ls observed in Fig. |1 
* 


Quasi-periodic Search and Analysis QPOs have been observed previously 
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during three giant flaring episodes from SGRs 0526-66, 1806-20 and 1900+14 
at frequencies ranging from 18-625 Hz [56] [57] [58] 59]. Of these, the QPO 
detections for SGRs 1806-20 and 19004-14 were restricted to the oscillating tails 
well after the conclusion of the initial spike. These periodicities have been 
interpreted as signatures of torsional sub-surface seismic oscillations triggered 
by the cataclysmic rupturing of the crust seeding the GFs. 

We searched for a spin frequency over a range of v — 0.02-50 Hz using both an 
unbinned and a logarithmically binned periodogram, in particular, we searched 
for signals with at least p « 0.01 (corrected for the number of frequencies and 
segments). We did not detect any signal that could be associated with stellar 
rotation. 

We searched for QPOs in the Nal GBM data for GRB 200415A using the 
same methodology as that used to search for QPOs in the GBM data for SGR 
J0501--4516[60], over a range of v — 40 — 4000 Hz. Due to telemetry packets 
dropping shortly after the peak, we searched three segments independently: 
the initial spike, Ts = [-5,5] ms, the fall time, T; = [7,160] ms, and a long 
segment of 200s after the initial spike to look for potential neutron star pulsation 
nominally in the range 0.05-1 Hz in the putative sub-background emission. We 
find no credible signal in the initial spike or the long segment. In the burst 
decay, we find a potential broad QPO candidate with moderate significance at 
v ~ 180 Hz. As the signal is broad, we use three different strategies to establish 
its significance. First we determine the trial-corrected p-value (probability) for a 
logarithmically-rebinned periodogram under the assumption of pure white noise 
is p = 8.3 x 107^. However, the candidate's frequency places it at the edge of 
the frequency regime where stochastic variability in the form of red noise is 
important. To test whether the QPO can be explained by a combination of the 


overall decay of the light curve and the associated window function, we fit the 
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light curve of this segment with an exponential function, and draw 1000 sets of 
parameters from a multi-variate Gaussian using the best-fit solution and inverse 
Hessian derived from the optimization. We use these parameters to generate 
simulated light curves from the exponential function by adding white noise and 
generating periodograms in the same way we did for the data. The p-value for a 
QPO at that frequency under the assumption of an exponentially decaying light 
curve is p — 0.023, indicating that while less decisive, a decaying exponential 
cannot alone explain the observed excess in power at 180 Hz. We then use the 
methodology from Huppenkothen et al.[60] and model the periodogram with a 
power law (low-frequency red noise) and a constant (white noise), then search 
for outliers in the residuals. Using this method, we find a posterior-predictive p- 
value of p — 0.016, indicating that pure red noise cannot easily explain the excess 
power in the observed data. Finally, as the QPO is fairly broad, we model the 
periodogram with a combination of a power law, a Lorentzian and a constant, 
then compute the likelihood ratio between this model and the red noise-only 
model. Using the latter, we generate red noise via Markov-Chain Monte Carlo 
simulations using emcee and derive a posterior predictive p-value for the 


likelihood ratio, p = 0.031. The parameters are well-constrained and uni-modal. 


We find a posterior centroid frequency for the QPO of v = 183.2+19.9 Hz, with 


a width of y = 46.7 + 29.7 Hz. We used the posterior distribution of QPO 
parameters to derive the fractional root mean squared (RMS) amplitude, and 


find Aga, = 0.39 + 0.017. Overall, while the significance is not consistently 


above a 3c threshold, the significance remains consistently small in all attempts 
of quantification. 

We similarly searched for the same signal in data from the GBM BGO detec- 
tors, as well as BAT data and find no counterpart in either, which may be related 


to the lower count rates in the BAT, as well as a potential energy-dependence of 
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the candidate. Using the aforementioned methods, we performed independent 
searches (both in the GBM and BAT data) at and around the QPO frequencies 
reported by ASIM[I6], but did not find any QPO candidates. 


Radio Search Four observations of the Sculptor Galaxy were taken with the 
Karl G. Jansky Very Large Array (VLA) [20] under the project VLA/20A-438, 
4.3 to 51.2 days post burst with the full 1-2 GHz L-Band. 3C48 was used for 
the band-pass and flux calibration and J0025-2602 was used for the complex 
gain calibration. 

The flagging, calibration, and imaging were done using the Common Astron- 
omy Software Application (CASA) [62]. Automated flagging was done using the 
tfcrop and rflag routines. After flagging and calibration, the spectral windows 
that were free of radio frequency interference were selected for imaging and 
four rounds of self-calibration. The baselines shorter than 2.5 kilo-wavelengths 
were excluded from the imaging process due to the bright diffuse emission of 
the Sculptor galaxy. The resulting images have a FoV diameter of ~ 1/3 of a 
degree. The dates of each 11 minute observation, along with its accompanying 
RMS noise are listed in Table 2] We note that due to the bright emission from 
the center of the galaxy, there are some artifacts from the VLA arms affecting 
the RMS noise in parts of the image by a factor of ~ 1.5. Overall, the noise 
in these images is high due to the bright point sources in the image and the 
removal of short baselines due to the bright extended emission of the galaxy. A 
direct comparison between images, as well as image subtraction using the CASA 


immath task, shows no significant variable or transient emission in the field. 


Data Availability Gamma-ray data from CGRO-BATSE, Swift-BAT and 
Fermi-GBM can all be found in public repositories. Catalogs of these data 
are provided as citations in both the main text and the Methods. The raw VLA 


data is publicly available. The calibrated VLA data and images are accessible 
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to the reader per request. 


Code Availability Standard software packages such as 'rmfit for GBM and 
‘XSPEC for other instruments, are available on-line. The codes used to deter- 
mine the significance of the BGO photons, to construct the BAT TTE DRMs 
and to determine the rise-time, are available per reasonable request. The VLA 
data was analyzed with publicly available software (CASA). The entire pro- 
cedure for detecting and quantifying the QPOs implemented in this paper, is 
publicly available in ‘Stingray’. The algorithm used to determine the pulse pile- 
up of the GBM data is available in Vandiver et al., (2013) [40]. The SwiMM code 
is not publicly available. However, response functions can be used to reproduce 


the spectral results in this study and are available upon reasonable request. 
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ED Fig. 1. The duration of GRB 2004154. The Too (green) and T5o (purple) 
durations, calculated using the Swift-BAT data in counts-space. The errors are 
at the l-o confidence level. 
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ED Fig. 2. Spectra and fitted models in three time intervals for GRB 

200415A. vF, spectra and Comptonized fitting residuals for Interval 1 (Left), 

Interval 3 (Center) and Interval 4 (Right) of GRB 2004154. The three spectra 
are devoid of any instrumental effects attributed to band-width saturation, 

with the fit parameters listed in Table[I] 'These figures show the robustness of 
the fits to the data (l-o confidence) which are used in the main text and in 

Fig. Ep, and are a direct result of the unrivalled temporal and spectral quality 

of the GBM data. 
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ED Fig. 3. Energetic photons from GRB 2004154. Individual 
Time-Tagged Events of GBM BGO detector 0 (black dots). The blue rectangle 
indicates energies from 2.5 to 3.5 MeV in intervals (2) and (3), while the red 
rectangle shows energies from 3.5 to 10 MeV. We conclude that the highest 
photon energy unambiguously associated to GRB 200415A is 3 MeV. 
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Date AThurst | lo RMS noise 
(UTC) (days) | (mJy/beam) 
April 19, 16:12:36 | 4.31 0.28 
April 25, 16:26:53 10.3 0.43 
May 7, 15:04:53 42.3 0.40 
June 5, 12:43:29 51.2 0.29 


Table 2: 1-2 GHz VLA radio observations 
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